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Abstract 

The Schrodinger equation discribing the local potential model of a strong tt + tt~ -interaction was studied. The 
influence of the strong ^^"-interaction of the behaviour of pionium nS-state wave functions at small distances is 
studied both analytically (perturbatively) and numerically. It is shown that in the whole the accounting of strong 
interaction results in multiplying pure Coulomb pionium wave functions by some function is practically independent 
of the value of the principal quantum number n. Due to this reason, the n-independence of probability of 7r + 7r~-atom 
production in nS-state remains the same as in the case of a pure Coulomb 7r + 7r~-interaction. 
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1. Introduction 



The method for the measurement of the pio- 
nium (7r + 7r~-atom) lifetime proposed in ref. jy, 
is essentially based on the assumption that the in- 
dependence of probability w n on creation of 7r + 7r~- 
atom in raS-state is well known or, at least, may be 
calculated with a high degree accuracy. The first 
consideration of this problem has been done in L. 
Nemenov's paper j!| , where the following relation 

w n ~ n~ 3 (1) 

has been derived from a more general result of the 
author § 



M(r 



M(r)^ n (r)d 3 r 



(2tt) s 



M(p)e- l P p d 3 p, 



(2) 



(3) 



where M(p) is the amplitude of production of free 
7r + 7r~ -pairs with relative momentum p in hadron - 
hadron or hadron - nucleus collisions and ip n (r) is 
the wave function of nS-st&tc of pionium. 
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In line with accounting the short range nature 
of amplitude M(r) in his original derivation of (1) 
from (2) L. Ncmcnov also has made an assump- 
tion that the pure Coulomb wave functions describe 
quite well a distribution of pions in the pionium 
not only at large distances, but at small ones also. 
However, as it has been shown recently by E. Ku- 
raev this assumption is unjustified due to some 
noticeable influence of strong a 7r + 7r~-interaction 
on the behaviour of the pionium wave functions in 
the nearest of origin. Due to this reason a more 
careful analysis of this problem is needed. 

Below we shall represent some preliminary re- 
sults of the analysis based on the local potential 
model of the strong 7r + 7r _ -interaction. In this 
model the "reduced" pionium wave functions 



(4) 



$n(r) = V47rr^„(r) , / $„(r) \ dr = 1 



obey the Shrodinger equation 

K(r)+m[U c (r) + U s (r)]$ n (r) = ms n $ n (r) , (5) 

where m is the pion mass, e n — binding energy, 
U c = a/r, U s — Coulomb and strong potentials, 
respectively. 
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2. The scheme and results of perturbative 
analysis 

First of all, let us apply the methods of the per- 
turbation theory to the Schrodinger equation (5), 
treating the strong interaction potential U s as per- 
turbation, in order to obtain some qualitative esti- 
mations. Putting 



$„(r)«^°>(r) + $«(r), 
where 



F - P (o)+ F (i) 



r (0) 
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W= / U s (r)^ 2 (r)dr; 



4°)"(r)+m[f/ c (r)-4 )]$( t °)(r) = 0, 



S«"(r) + m[C/ c (r) £ (°)]$W(r) 



(G) 
(7) 

(8) 

(9) 



m[ e W-I7.(r)]*f(r) 



(10) 



and applying the general methods of solving the 
linear inhomogeneous equations one can ob- 
tain 



where 



<^ 0) (n)| 



(ii) 



(12) 



x [- £ W + f/ s (r 2 )]dr 2 



and 



c„= / dr|*<°>(n)| \n{r). 
o 

If we define the ratio 



Rn{r) 



(c)/ 



(13) 



(14) 



as a measure of the influence of strong interactions 
on the values of the pionium wave functions, then 
in the first order of the perturbation theory 



Rn(0) PS 1 + C n . 



(15) 



With explicit expressions for the pure Coulomb 
wave functions $1°^ we have proceeded in the cal- 
culation of c n with n = 1, 2, 3. 



The result looks like as follows 

An) 



x 



71=1 

fc / (n) • 



™ =1 

OC ^ 



where <i^™' 1 ~ 1, r^™' ~ w 400 /m. 
Taking into account the relation 



mU s (r)r 2 dr ps a w 0.15 /m 



(16) 



(17) 



where a is the 7r + 7r~ scattering length, it is easy 
to see that the n-dependent contributions to c n are 
numerically small (of order 10 -3 ) and may be ne- 
glected. 

Putting U s = g/r exp(— 6r) with the values of 
parameters f| g ps 3, b = m p ps 3.8 fm^ 1 , that 
corresponds to applying the p-exchange model for 
a describtion of the strong 7r + 7r~-interaction, one 
can obtain for n-independent part of c n (see also 



1) 



mU s (r)rdr 



gm 



0.55. 



(18) 



These estimations show that the strong ir + ir~- 
interaction can change noticeably the value of the 
pionium wave functions at origin and this effect 
cannot be ignored at evaluating the probabilities 
w n (2). On the other hand, the large correction to 
the values of the wave functions, obtained in the 
first order of perturbation theory, means that the 
higher order corrections are not small and must be 
taken into account. Their calculation does not look 
a simple problem. 

3. Numerical investigation 



Due to this reason we have applied numerical 
methods for an accurate investigation of the pio- 
nium wave functions behaviour at small distances, 
using the following equation with boundary condi- 
tions and normalization: 

fx 

dp 2 



$ 1 = 7 4 + c/(p)xW-ex(p) = o, 



$ 2 = x(o) = o, $ 3 = xM = o. 



(19) 
(20) 



$4 = / X (P)dp 
a 



1 = 0, where 



(21) 
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Fig. 1. A result of numerical calculation of Coulomb 
pionium eigenfunction Xi^\p)- 
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So we have an operator equation 



(22) 



$(z) = : 



where 



(23) 



$(z) = < 



$2(2) 
$4(2) 



2 = (x(p), e) • 



(24) 



In order to calculate the normalized pionium 
wave functions, we have used an improved version 
of code Jt| , based on applying the continuous anal- 
ogy of Newton method, developed in ref. [||. In 
this way we change the primary equation by the 
following differential equation with parameter t 



^* (*(*)) = -*(*(*)) , 



z(0) = Zo, < t < 00 . 



(25) 



(26) 



When t — > 00, the solution of this equation tends 
to the solution of primary equation (23) z(t) — > z. 

For numerical calculations the grid on t was used 
{U, i = 1,2,... ; t ; £t+i -U = r}. For every tj 
we obtained the solution dz(t = U)/dt to equation 
(25) 



^^'(.(i,))- 1 $(*&)), 



(27) 



where 0' is Frechet derivative of function <j>'~ 
is an inverse operator. 

For approximation of dz(ti)/dt we use a finite 
difference 



dz(U) 

dt 
Thus, 



Az,; 



[*(ti+i) - z(U)] 



z{U+l) = TiAzi + Z(ti) 



(28) 



(29) 



If z(ti) is known, then as a result of this it- 
erative procedure with parameter n we obtained 
a sequence of the approximate values of the so- 
lution to equation (25) {zi} — > z. As a start- 
ing approximation zq we choose an analytical solu- 
tion of the corresponding Coulomb problem. The 
iterative procedure was continued until residual 
||$(z(i))||<10- 5 . 

The problem (27) was solved numerically in the 
interval (a, b), where a = 10~ 8 , & is a sufficiently 
large number, depending on the principal quantum 
number n. The solution at the point b was ad- 
justed logarithmically to the appropriate Coulomb 
eigenfunction 

X'(b) x {c) '(b) 



X(b) X^Kb) 



(30) 



To calculate, we used a nonuniform grid on p as 
follows: in the interval (a, 0.02) — 301 nodes, in the 
interval (0.02, 2) — 301 nodes and in the interval 
(2, b) — 200 nodes. 

The accuracy of the method was tested on 
Coulomb problem. Fig. 1 shows the numerically 
obtained solution of Coulomb problem for n = 1 
for potential U c = 2/ p. Fig. 2a gives the corre- 
sponding value of the wave function ip^ (p) . The 
same functions for n = 2, 3, 4 are given in Figs. 
2b-2d. The difference between the numerical and 



analytical values of function ip. 
is represented in Figs. 3a-3d 
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Fig. 2. Results of numerical calculations of Coulomb 
pionium nS-state wave functions vi c '(p) for n = 1 (a), 
n = 2 (b), n = 3 (c), n = 4 (d). 
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Fig. 3. Difference between numerical and analytical val- 
ues of functions ipn\p) for n = 1 (a), n = 2 (b), n = 3 
(c),n = 4(d). 
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(a) 



Fig. 4. Behaviour at small distances on the pionium 
wave function tpi(p). 

The input of parameters of the code has been cho- 
sen in this way to guarante an absolute accuracy of 
the calculations higher than 1CP 4 . To check up this 
accuracy, we have compared the numerical solution 
of Shrodinger equation with the pure Coulomb in- 
teraction with the analytical one. 

At the next step we solved our problem for the 
sum of strong and Coulomb potentials as follows: 



U(p) = -(l-e- b P) + -e- b P. 



where 



(31) 



fe=!^ s ,1.5.10 3 , o = 
ap 



2a 



8 • 10 2 , (32) 



Q 2 

v pirn 
47T 



3. 



The results of the numerical calculations of the 
wave function ip(p) for n = 1 are represented in 
Fig. 4. At a long interval of p we cannot see a 
difference between it and the corresponding wave 
function of Coulomb problem, but when p is com- 
patible with Fermi radius, then this difference is 
significant. This fact is shown for n = 1 in Fig. 5a. 
The difference also takes place for n = 2, 3, 4 (see 
Figs. 5b-5d). 
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Fig. 5. Influence of strong pion-pion interaction on the 
pionium nS-state wave functions tpn(p) at small dis- 
tances : n = 1 (a), n = 2 (b), n = 3 (c), n = 4 (d). 



The most important conclusion from our numer- 
ical calculations is as follows: the functions 



MP) = (33) 
V-n (P) 

are practically independent of n (n — 1, 2, 3, 4). 
This is illustrated by Fig. 6. 

Fig. 7a shows the differences between functions 
R n (p) for the considered n. One can see from the 
plots that the differences between functions R n for 
n = 1 — 4 are of the order less or equal 1CP 3 . Thus, 
we can replace functions R n (p), {n = 1 — 4) with 
accuracy 10~ 3 with R±(p) 

R n (p) = Riip) + O(10- 3 ) . (34) 
The same results were obtained for potential 

U(p) = - (1 - e- bp ) + a e- bp , where (35) 
P 

b= Utt rj 1.5. io 3 , 5=i a -6« 6-10 5 . (36) 
a/i 2 



As one can see from Fig. 7b, the estimation (34) is 
valid for potential (^5|). 

So, we can make the following conclusions from 
our calculations: for the considered potentials ( |3l| ) 
and (|35|) for the principal quantum numbers n 
1 . 2, 3, 4 the following estimation is valid 

^ - + o flO- 3 ) (37) 

i.e. with accuracy 10 -3 function ifi n (p) can be rep- 
resented as 

Mp)*u(p)^n C) (p)i (38) 

where a(p) is independent of n, p is compatible 
with Bohr radius. 



4. Conclusion 



Fig. 6. The ratio Rn(p) as a measure of the influence 
of strong interactions on the values of the pionium wave 
functions: n = 1 (a), n — 2 (b), n = 3 (c), n — 4 (d). 



Thus, the results of the numerical solution to 
equation (5) with Yukawa-type strong potential 
have confirmed the main conclusions of the per- 
turbative consideration, namely: the ratios 

Rn(r) = ^-^R n (p) (39) 
4>n (r) 

being numerically large (see Fig. 6) in the region 
r < r s ^ 1 fm and essential for the problem un- 
der consideration (n-dependence of values w n ), are 
practically n-independent (their n-independence is 
illustrated by Fig. 7). 



Fig. 7. The accuracy of n-independence of strong renor- 
malization factors for potentials (a) U(p) and (b) U(p). 



This means that with a high degree accuracy one 
can substitute 

i> n {r) =R{r)^{r) (40) 

in eq. (2) and, replacing 

M(f) M{r) = M(r)R(r), 



obtain 




Therefore, we can conclude that strong inter- 
action corrections to the pionium wave functions 
at small distances, being sufficiently large, do not 
change the n _3 -law (1) primarily derived in paper 
|!| assuming that the pionium wave functions are 
pure Coulomb. 
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